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Abstract 

We have studied the electronic structure of 3d transition-metal-Mgf^ alloys, Mgo.97TMo.03B2, 
(TM = Sc, Ti, V, Cr Mn, Fe, Co, Ni, Cu, Zn) using Korringa-Kohn-Rostoker coherent-potential 
approximation (KKR-CPA) method in the atomic-sphere approximation (ASA). For unpolarized 
calculations, our results for Mgo.97TMo.03B2 alloys are similar to that of 3d impurities in other s 
and s — p metals. In particular, the local densities of states (DOS) associated with the 3d impurities 
are similar to our earlier work on 3d impurities in bulk Al [P. P. Singh, Phys. Rev. B 43, 3975 (1991), 
P. P. Singh, J. Phys.: Condens. Matter 3, 3285 (1991)]. For spin-polarized calculations, we find only 
the alloys of V, Cr, Mn, Fe and Co with MgBi to be magnetic of all the 3d elements. We also find 
that Cr and Mn in MgF$2 have a relatively large local magnetic moment of 2.43 ub and 2.87 /j,b, 
respectively. We have used the unpolarized, self-consistent potentials of Mgo.97TMo.03B2 alloys, 
obtained within the coherent-potential approximation, to calculate the electron-phonon coupling 
constant A using the Gaspari-Georffy formalism and the superconducting transition temperature T c 
using the Allen-Dynes equation. We find that the calculated T c is the lowest for Mgo.97Cro.03B2 and 
the highest for Mgo.97Zno.03B2, in qualitative agreement with experiment. The calculated trend in 
variation of T c from Mn to Zn is also similar to the available experimental data. Our analysis of 
the variation in T c , in terms of the DOS and the spectral function along Tto A direction, shows the 
variation to be an interplay between the total DOS at the Fermi energy and the creation/removal 
of states along r to A direction [P. P. Singh, cond-mat /Q2Q 1093 1 . 



PACS numbers: 74.25.Jb, 74.70.Ad 
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I. Introduction 



The nature of interaction responsible for superconductivity in MgB 2 [|, |], £|, [5|, [7], 



10, 11, 12, 13, 14, O, 16 suggests a gradual decrease in the superconducting transition 



temperature, T c , upon addition of impurities [17], [18|, [T^, ^(J with increasing electron/ atom 
ratio. A systematic increase in the number of available electrons is expected to fill up the a 
holes in MgB 2 , which are coupled strongly to the in-plane bond-stretching mode of B, and 
thereby reduce the strength of the electron-phonon coupling resulting in a decrease in T c . 
The observed variation in T c of alloys of Al with MgB 2 and its understanding in terms of a 
gradual filling up of a holes provides a good example of this picture pij, p!fl . 

Since the electron /atom ratio increases as one goes from Sc to Zn, based upon the above 
argument one expects the T c to systematically decrease as the different elements from the 3d 
row are added to MgB 2 . However, the experimentally observed changes in T c of 3d transition- 
metal-M(7i?2 alloys |20j do not seem to follow the expected trend of a systematic decrease 
in T c as one goes from Sc to Zn. For example, a 3% Mn-doped MgB 2 (Mg 97 Mn 03 B 2 ) 
has a T c of only 33.1 K, while MgB 2 doped similarly with Fe, Co and Ni show T c 's of 
37.8, 35.7 and 37.8 K, respectively p0[| . Surprisingly, Zn— doped MgB 2 shows the highest 
T c (= 38.4 K) of all the transition-metal-Mg£> 2 alloys investigated so far (19|, ^0~| . The 
observed changes in T c of transition-metal-doped MgB 2 alloys, especially as one goes from 
Fe to Zn, are unexpected even after making allowances for magnetic effects [^3|, [24|| . 

To understand the changes in the electronic structure and the superconducting proper- 
ties of MgB 2 alloys upon addition of 3d transition-metal impurities, we have carried out 
ab initio studies of Mgo.g 7 TM .o 3 B 2 (TM = Sc, Ti, V, Cr Mn, Fe, Co, Ni, Cu, Zn) alloys 
using density-functional-based methods. We have used Korringa-Kohn-Rostoker coherent- 
potential approximation method [^] in the atomic-sphere approximation (KKR-ASA CPA) 
|26H for taking into account the effects of disorder, Gaspari-Gyorffy formalism |27|] for cal- 



culating the electron-phonon coupling constant A, and Allen-Dynes equation [23], |24| for 
calculating T c in Mg 0i g 7 TM om B 2 alloys. For understanding the variation in T c as one goes 
from Sc to Zn, we have analyzed our results in terms of the changes in the spectral function 
25j along r to A [21], ^8[ and the densities of states (DOS), in particular, the changes in the 



B p and the transition-metal d contributions to the total DOS. The changes in the magnetic 
alloys are described in terms of spin-resolved densities of states and the local magnetic mo- 
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ments. Before we describe our results, we outline some of the computational details of our 
calculation. 



II. Computational Details 

The charge self-consistent electronic structure of unpolarized as well as spin-polarized 
Mg 097 TM 03 B 2 alloys has been calculated using the KKR-ASA CPA method. We have 
used the CPA successfully to describe the electronic structure of Al— doped MgB 2 alloys 



2Tfl. We parametrized the exchange-correlation potential as suggested by Perdew-Wang 
within the generalized gradient approximation The Brillouin zone (BZ) integration 

was carried out using 1215 k-points in the irreducible part of the BZ. For DOS calculations, 
we added a small imaginary component of 1 mRy to the energy and used 4900 k-points in 
the irreducible part of the BZ. The lattice constants for Mgo, 97 TM om B 2 were fixed at the 
MgB 2 values. The Wigner-Seitz radius for Mg was slightly larger than that of B, while the 
Wigner-Seitz radii of the impurities were equal to their bulk value as given in Ref. |30]]. The 
maximum I used was l max = 3. 

As indicated above, the electron-phonon coupling constant A was calculated us- 
ing Gaspari-Gyorffy formalism with the charge self-consistent potentials of unpolarized 
Mgo.g?TMomB2 alloys obtained with the KKR-ASA CPA method. Subsequently, the 
variation of T c was calculated using Allen-Dynes equation. The average value of phonon 
frequency^ for MgB 2 was taken from Ref. |§ and //* = 0.09. 



III. Results and Discussion 



In this section we present the results of our self-consistent electronic structure calculations 
for Mg 097 TM 03 B 2 alloys. We first describe the results of the unpolarized calculations in 
terms of the total and the sub-lattice resolved DOS, including the local DOS due to the 
3d transition-metal impurity. From our spin-polarized calculations, described next, we find 
only the alloys of V, Cr, Mn, Fe and Co with MgB 2 to be magnetic of all the 3d elements. 
These results are discussed using the spin-resolved DOS and the local moments. Finally, 
the calculated variation in T c for Mg 097 TM , Q3 B 2 alloys is compared with the available 
experimental data, and analyzed in terms of the DOS and the spectral function along T to 
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.4. 



A. Unpolarized Total and Local Densities of States 

In Fig. 1 we show the calculated densities of states of Mg Q ,g 7 TM om B2 alloys from Sc 
to Zn, including the total DOS, N(E), and the concentration- weighted, sub-lattice re- 
solved DOS for the two inequivalent sub-lattices (N sub (E)) in MgB 2 . The two inequiv- 
alent sub-lattices consist of the Mg and the transition-metal on one sub-lattice (N{ ub (E) = 
0.97N s ^ h g (E) + 0.03N^(E)) and B on the other two sub-lattices (iV|" b (£) = 2N s B ub (E)). As 
we go from Sc to Zn, the increase in the electron/atom ratio moves the Fermi energy, Ep, as 
well as increases the hybridization of the host s — p bands with the impurity d— band. The 
change in Ep has important consequences for the superconducting properties of these alloys, 
as discussed later. The inward movement of the d-b&nd as we go from Sc to Zn is clearly 
manifested in the total DOS as well as the concentration-weighted DOS of Mg sub-lattice. 
However, we find that the DOS at the B sub-lattice remains largely unaffected due to the 
transition-metal impurity, although the movement of Ep changes the B contribution to the 
total DOS at Ep. 

To be able to examine the changes in the DOS of the transition-metal impurities and to 
show the movement of the d— band as we go from Sc to Zn, we show in Fig. 2 the DOS 
of the transition metal impurity N^(E) on the Mg sub-lattice in Mg ,g 7 TM om B2 alloys. 
We find that the d— level crosses Ep between Cr and Mn, resulting in a very high density 
of states at Ep for these alloys. The relatively high DOS for N^ b (E F ) and N s ^ b n (E F ) points 
to the possibility of local magnetic moment formation at the impurity sites. By the time we 
come to Zn impurity, the d— level is well inside Ep and the impurity contribution to the total 
DOS has reduced significantly. We also find that the lowering of d— level with increasing 
d— electrons is accompanied by the narrowing of the d— band. For example, in the case of 
Zn impurity the width of the d— band is the narrowest, indicating almost an atomic-like 
behavior of these d— electrons. It is interesting to note that these results are qualitatively 



similar to our results on 3d transition-metal impurities in Al [31 , |32 
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FIG. 1: The calculated total density of states (solid line) for Mgo.97TMo.03B2 alloys. The contri- 
butions to the total DOS from the Mg sub-lattice (dashed line) and the B sub-lattice (dotted line) 
are also shown. Note that the Mg sub-lattice contains both Mg and the transition-metal (TM) 
impurity atoms. The vertical dashed line indicates the Fermi energy. 
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FIG. 2: The calculated total local density of states (solid line) for transition-metal impurity in 
Mgo.97TMo.03B2 alloys. The vertical dashed line indicates the Fermi energy. 
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FIG. 3: The calculated total density of states (solid line) for M g^^T Mq^B 2 alloys with majority 
(upper panel) and minority spins (lower panel). The majority and minority spins contributions to 
the total DOS from the Mg sub-lattice (dashed line) and the B sub-lattice (dotted line) are shown 
in their respective panels. The vertical dashed line indicates the Fermi energy. 

B. Spin-Polarized Total and Local Densities of States 

The results of our spin-polarized calculations for Mg OQ7 TM om B 2 alloys are shown in 
Figs. 3-5. Our calculations show that the alloys of V, Cr, Mn, Fe and Co in MgB 2 are 
exchange-split, as can be seen in Fig. 3 where we have plotted the spin-polarized total and 
the concentration-weighted, sub-lattice resolved DOS for V, Cr, Mn, Fe and Co in MgB 2 . 
Once again, we find that the B sub-lattice remains unaffected due to magnetic moment 
formation at the impurity site on the Mg sub-lattice. To clearly show the changes in the 
majority and the minority spin DOS of the Mg sub-lattice due to exchange-splitting in 
the DOS of the impurity atom, we show in Fig. 4 the spin-resolved total DOS due to the 
impurity atom. As expected from the unpolarized calculations, the exchange splitting is the 
largest for Cr and Mn alloys. The exchange-splitting in the DOS due to the impurity atom 
leads to local moment formation at these sites with Cr and Mn having local moments of 
2.43 /is and 2.87 //b, respectively. The local moment on Co is very small (0.01 hb)- Note 
that the total DOS of the impurity is dominated by the d contribution and the local moments 
arise almost entirely due to the impurity d electrons. 



6 




1 -0.5 -1 -0.5 -1 -0.5 -1 -0.5 -1 -0.5 
E(Ry) 



FIG. 4: The calculated total local density of states (solid line) for transition-metal impurity in 
Mgo.g7TMo.03B2 alloys with majority (upper panel) and minority spins (lower panel). The vertical 
dashed line indicates the Fermi energy. 



It must be pointed out that the calculation of local moments in alloys is sensitive to 
the volume that one associates with the impurity atom. Given the availability of electrons 
around Mg sub-lattice in MgB 2 alloys, a judicious choice of the impurity-atom volume is 
essential for a reliable description of the local moments. In our calculations, we have chosen 
the impurity atomic volume to be equal to the observed bulk value as given in Ref. [pCjll. 
However, it is clear that Cr and Mn will show large local moments in MgB 2 . 



C. Superconducting Transition Temperature 

We have studied the electronic structure of Mg 97 TM 03 B 2 alloys as a prelude to un- 
derstanding the superconducting properties, especially the superconducting transition tem- 
perature T c , of these alloys. In Fig. 6 we show our calculated T c as well as the observed 
T c [2(J for Mgo.97TMo.03B2 alloys. The calculated variation in T c across the 3d row is 



similar to the one observed experimentally for Mn to Zn in MgB 2 (T^, [2IJ. Note that our 
calculated T c for MgB 2 is equal to ~ 30 K, which is consistent with the results of other 
works 0, P, |l(J with similar approximations. 

The total DOS and the spectral function along r to A at E F are expected to play an im- 
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MgB 2 Sc Ti V Cr Mn Fe Co Ni Cu Zn 



FIG. 5: The calculated local magnetic moment (filled circle) of the transition-metal impurity in 
Mg , 97 TM om B 2 alloys. 

portant role in deciding the T c of Mg 97 TM Q 3 B 2 alloys. Thus, in the following we examine 
the changes in these two quantities in Mg 097 TM 03 B 2 alloys as we go from Sc to Zn. In Fig. 
7 we show the total DOS, N(Ep), the Bp contribution to N^ b (Ep) and the impurity d con- 
tribution to Ntm(Ef) at the Fermi energy. Our calculations show Mg , 097 Croxi 3 B 2 to have 
the lowest T c of all the 3d alloys, which coincides with the highest N(Ep) of 11.93 stj (Ry cell) 
as well as the highest d contribution to the Nq^(Ep) of 66.15 st/(Ry atom), as can be seen 
in Fig. 7. In contrast, the Bp contribution is not enhanced resulting in the lowest T c for 
Mg Q97 Cr 0Q3 B 2 within our approach. Similarly, the T c for Mg , 97 Mn ,o 3 B 2 is also small in 
comparison with other alloys. Here, we like to point out that the exchange-splitting will 
reduce Nj^(Ep), leading to a smaller N(Ep) in the case of magnetic alloys. Thus, within 
our approach, it would have led to an increase in T c . However, the inclusion of magnetism 
with its pair-breaking effects p3l B|| will further reduce the T c . The gradual increase in T c 



from Mn to Zn is due to enhanced Bp contribution to Ng" b (E F ) as well as a substantial 
decrease in the impurity d contribution to N^(Ep) as shown in Fig. 7. 

To further clarify the reasons for the changes in T c in Mg .g 7 TM om B 2 alloys, as we go 
from Sc to Zn, we show in Fig. 8(a)-(d) the spectral function, A(k,Ep), calculated along T 
to A direction in Mg . 97 TM om B 2 alloys. For comparison we have also plotted A(k,Ep) for 
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MgB Sc Ti V Cr Mn Fe Co Ni Cu Zn 



FIG. 6: The calculated (filled circle) and the experimental (filled square) variation of T c in 
Mg , 97 TM .o 3 B 2 alloys. 



MgB 2 in Figs. 8(a) and (d). As described in Ref. pjj , the additional states created along 
r to A direction reduce the coupling between the hole-like cylindrical Fermi sheets with the 
phonons which, in turn, reduces the T c . Thus an alloy with A(k,E F ) similar to that of MgB 2 
will show a T c close to that of MgB 2 . Carrying out such a comparison in Fig. 8, we find 
that Mgo.g 7 Zn om B 2 will have a T c close to that of MgB 2 , while Mgo.$ 7 Vo,o3B 2 will show the 
lowest T c of them all. The fact that our calculation shows Mg , 97 Cro,o3B 2 to have a lower 
T c than Mg 097 V om B 2 is due to the relatively high impurity DOS at E F . In addition, in the 
case of Cr and Mn impurity in MgB 2 , the magnetic effects will further reduce the T c p4|. 



These results reinforce the observation made in Ref. [gT| that the way T c changes in MgB 2 
upon alloying depends dramatically on the location of the added/modified k-resolved states 
on the Fermi surface. 



IV. Summary 

We have studied the electronic structure of 3d transition-metal-Mg£> 2 alloys using 
Korringa-Kohn-Rostoker coherent-potential approximation method in the atomic-sphere ap- 
proximation. Our results for the unpolarized calculations are similar to that of 3d impurities 
in other s and s — p metals. From spin-polarized calculations we find that only alloys of 
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MgB 2 Sc Ti V Cr Mn Fe Co Ni Cu Zn 



FIG. 7: The calculated total density of states at the Fermi energy (filled circle, upper panel), Bp 
contribution (filled circle, middle panel) and the transition-metal (TM) d contribution (filled circle, 
lower panel) to the total DOS in M go.^T 'M .03-B 2 alloys. 




FIG. 8: The calculated spectral function along T to A, evaluated at the Fermi energy in 
Mgo.97TMo.03B2 alloys. In figures (a)-(d) the solid, dashed and dotted lines correspond to (a) 
Sc, Ti and MgB 2 , (b) V, Cr and Mn, (c) Fe, Co and Ni, and (d) Cu, Zn and MgB 2 alloys, 
respectively. Note the change in the vertical scale in (6). 
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V, Cr, Mn, Fe and Co in MgB 2 are magnetic, with Cr and Mn having the largest lo- 
cal magnetic moment of 2.43 hb and 2.87 /xb, respectively. We have used the unpolarized, 
self-consistent potentials of Mgo.g7TMo.Q3B2 alloys, obtained within the coherent-potential 
approximation, to calculate the electron-phonon coupling constant A using the Gaspari- 
Georffy formalism. Then, with the help of Allen-Dynes equation we have calculated the 
superconducting transition temperature T c of these alloys. We find that the calculated T c is 
the lowest for Mgo.g7Cro.03B2 and the highest for Mgo.g7Zno.03B2, in qualitative agreement 
with experiment. The trend in variation of T c from Mn to Zn is also similar to the available 
experimental data. 
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